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ABSTRACT 

Turbulence structures in a wall-bounded shear layer during the 
bursting event detected by a conditional sampling technique are investi- 
gated using data obtained from large-eddy simulation of turbulent 
channel flow. Streamlines are constructed from the ensemble-averaged 
velocity field to illustrate the flow patterns associated with the 
bursting event. They exhibit the splatting motions during the sweep 
event and the existence of a pair of counterrotating streamwise vor- 
tices during the ejection process. 
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I. INTRODUCTION 


Organized structures in turbulent flows have been a subject of much 
investigation in the past two decades. Despite a great deal of effort 
to identify their characteristics, the lack of consensus on the detailed 
description of the structures still prevails. Comprehensive reviews on 
the subject are recently given by Cantwell 1 and Hussain. 2 

In this paper, detailed flow patterns associated with the bursting 
event in the near-wall region of turbulent boundary layers are investi- 
gated. The bursting phenomenon is composed of a sequence of quasi- 
cyclic events that occur in the wall region of turbulent flows. During 
this event, approximately 70% of total turbulence production occurs 
(see Kim et al„ 3 and Kline 4 ); hence, this is regarded as one of the 
most important processes in wall-bounded turbulent shear flows. An 
understanding of this phenomenon is essential to a better insight into 
the dynamics of the turbulent transport process. Because of its impor- 
tance, the bursting phenomenon has attracted a great deal of theoretical 
and experimental interest. Nevertheless, the exact description of this 
event, as well as its cause, is still in question. Kline 4 and Brodkey 5 
provide excellent perspective views based on their pioneering work 
accumulated in over a decade at Stanford and Ohio State universities, 
respectively. Recently, Smith 6 proposed a synthesized model derived 
from his flow-visualization work. Most of the previous work, however, 
fails to provide a clear picture of flow patterns associated with the 
bursting event because of the difficulty involved in obtaining the 
necessary information experimentally. Hussain 2 summed up the difficulty 
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as follows: "Flow visualization presents excessive information but 

very little hard data, and anemometer data give some hard data but 
very limited flow physics." Recently, with the advent of large and 
fast computers, computer simulations of turbulent flows have shown that 
these simulations can be used as alternative research tools (see Rogallo, 7 
Mom and Kim, 8 Kim, 9 and Moin, 10 for example). From these time- 
dependent three-dimensional simulations, one can obtain any desired 
information on the instantaneous three-dimensional velocity and pressure 
field. 

The study reported here is a continuation of the author’s previous 
work, 9 m which the conditional sampling technique developed by 
Blackwelder and Kaplan 11 was applied to computer-generated data obtained 
from a large-eddy simulation of turbulent channel flow. 8 Various 
ensemble-averaged quantities, such as velocity, pressure, and vorticity, 
were discussed in conjunction with the bursting phenomena detected by 
the conditional sampling technique. In the present work, more detailed 
flow patterns are described in terms with streamlines and contours con- 
structed from the three-dimensional ensemble-averaged velocity field. 

In the aforementioned work of Moin and Kim, a fully developed 
turbulent channel flow was simulated numerically at a Reynolds number of 
13,800, based on the centerline velocity and the channel half-width. 

The large-scale flow field was obtained by integrating the filtered, 
three-dimensional Navier-Stokes equations. The small-scale field motions 
were simulated through an eddy-viscosity model. The calculations were 
carried out with 516,096 grid points. The grid spacings m the 
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streamwise and spanwise directions were Ax ~ 62 and Az 2 15 in wall 
units. Nonuniform meshes were used in the normal direction, and the 
first mesh point away from the wall was y + « 1.78. The results were 
verified by comparing the statistical properties of the computed flow 
field with experimental results of turbulent channel flows. In addition, 
the turbulence structures in the vicinity of the wall-of-the-channel 
flow were found to be the same as the structures in the wall region of 
turbulent boundary layers. These data base obtained from the channel 
simulation were used in Kim 9 and the present work to study the flow 
structure associated with the bursting event in the near-wall region 
of turbulent boundary layers. 

II. CONDITIONAL SAMPLING 

To obtain an ensemble-averaged flow structure of the bursting 
event, the variable-interval time-averaging (VITA) technique developed 
by Blackwelder and Kaplan 11 is chosen as the conditional sampling 
process. In the present study (as in Kim 9 ), however, the conditional- 
averaging process is slightly modified to obtain a spatial structure 
rather than the temporal structure as in Blackwelder and Kaplan; the 
structures can be approximately related to each other by using 
Taylor’s hypothesis. It should be noted that most of the quantitative 
measurements in a laboratory experiment produce temporal structures , 
whereas most flow-visualization experiments generally yield spatial 
structures. A summary of the conditional sampling technique employed 
in the present work is given below, but the reader should refer to 
Blackwelder and Kaplan 11 and to Kim 9 for details. 
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The variable-interval space-averaging (VISA) is defined as 
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where L is the width of the spatial averaging. For all the data 
discussed in this paper, a threshold value, k = 1.2 and L = 8Ax, 
corresponding to 500 wall units (i.e., l"*" = Lu /v = 500), were used. 

T 

The detection criteria were applied at y + = 21. These procedures were 
applied to the data base generated by Moin and Kim, 8 and the results 
reported in this paper were obtained by ensemble-averaging over more 
than 3000 independent events. 


III. TURBULENCE STRUCTURE 

To illustrate the flow patterns associated with the bursting 
event, streamlines are constructed from the ensemble-averaged velocity 
field in several different planes. Figure 1 shows streamlines m the 
(x-y) plane at z = 0, where x, y, z correspond to the streamwise, 
normal-to-the-wall , and spanwise directions. These streamlines are 
constructed in a frame of reference moving at a speed of 0.83 U^ , 
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where U c is the mean velocity at the centerline of the channel. This 
speed was chosen as an assumed convection speed of large-scale structures 
(the correct convection speed of the large-scale structure is not 
known, but the qualitative picture does not change with different con- 
vection speeds). As a reference, the ensemble-averaged streamwise velocity 
at y + = 21, where the detection criteria were applied (x + = 0, y + = 21, 
z* = 0) , is given at the upper half of the figure. Note that the stream- 
wise extent of the figure is about 4000 wall units, and the vertical 
extent is from the wall to the centerline of the channel, which corresponds 
to 640 wall units. At this convected frame of reference, there exists 
a saddle point just upstream of the detection point exhibiting a similar 
flow pattern as in a turbulent mixing layer. 12 

Streamlines in the (y-z) planes, which are perpendicular to the 
mean flow direction, are shown in Figs. 2 and 3. The horizontal extent 
of these figures is about 600 wall units. Streamlines in Fig. 2 are at 
the upstream of the detection point, x + ~ -124, and show that high-speed 
fluid moves toward the wall (sweep event), thus creating an excess in 
the streamwise velocity, as indicated in the upper part of Fig. 1. 
Downstream of the detection point, x + = 124, Fig. 3 reveals that under- 
neath the sweep event a pair of counterrotating streamwise vortices is 
formed such that between them the low-speed fluid from the wall region 
is pumped away from the wall (ejection event) , creating a defect in 
the streamwise velocity. The corresponding velocity vectors are 
plotted in Figs. 4 and 5. Here, the velocity vectors are drawn parallel 
to the vectors formed by the spanwise and normal component of the 
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velocity. The purpose of these vector plots is to check the assumption 
one has to make in order to construct the streamlines in the (y-z) 
planes — that is, 3u/9x ~ 0. By comparing Figs. 2 and 3 with Figs. 4 
and 3, we notice no discernible difference in the flow patterns; 
this validates the approximation used in constructing Figs. 4 and 3. 

Contours of streamwise vorticity of the ensemble-ave raged velocity 
field m (y-z) planes are shown in Figs. 6 through 8. The sense of 
rotation of the vorticity is such that positive vorticity (solid lines) 
indicates clockwise rotation, and negative vorticity (dashed lines) 
indicates counterclockwise rotation. In Fig. 6, upstream of the 
detection point x + ~ -124, we notice that the sweeping motion of the 
high-speed fluid, depicted as a pair of counterrotating vortices in 
the figure, induces vorticity of opposite sign near the wall owing to 
the no-slip boundary condition. This viscosity-induced vorticity 
appears to be the source of the counterrotating streamwise vortices 
that appear downstream of the detection point, as shown m Figs. 7 
and 8. As the streamwise vortices are moving away from the wall 
because of the self-induced motion as well as because of the imposed 
pressure gradient (see Kim 9 for details) , they not only destroy the 
original vortices above, but also induce the opposite sign streamwise 
vorticity underneath them. Note also that the cores of the stream- 
wise vortices approach each other as they move away from the wall. 

In Fig. 7, the distance between the centers of the counterrotating 
vortices is about 120 wall units and the centers are located about 
23 wall units away from the wall; they are 75 wall units apart at 
about 45 wall units away from the wall m Fig. 8. (The scaling of 
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the large-scale structure was not captured properly in the simulation 
of Mom and Kim 8 because of the numerical resolution problem; therefore, 
any quantitative information must be taken with some care.) 

Contours of streamwise and normal velocities in the (x-y) plane 
are shown in Figs. 9 and 10. In these figures, solid lines represent 
positive values and dashed lines represent negative values. For the 
streamwise velocity, the mean value is subtracted from the ensemble- 
averaged velocity; hence, the positive and negative values are relative 
to the local mean value. These figures show that the high-speed fluid 
moving toward the wall is riding on top of the low-speed fluid in the 
wall region. Figures 11 and 12 show the contours of the streamwise 
velocity m the (x-z) plane at y + = 12 and (y-z) plane at 
x + = 186. These figures show that the low-speed fluid is completely 
surrounded by the high-speed fluid except at the downstream end, 
and it appears that the high-speed fluid causes the low-speed fluid to 
be squirted out. These figures are consistent with the description of 
Praturi and Brodkey 13 that the low-speed fluid is trapped between fingers 
of the high-speed fluid. 

Recently, Smith, 6 in an attempt to illustrate the existence of 
horseshoe vortices m the near-wall region of turbulent boundary layers , 
demonstrated that there exists a strong similarity between the laminar 
flow over a hemispherical protuberance and a turbulent flow during 
the bursting events. He proposed a model, in which the horseshoe vor- 
tices are formed m the low-speed streak regions by a process similar to 
that in laminar flow where the horseshoe vortices are formed when 
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vortex sheets are peeled off from the hemispherical surface. There is 
no evidence in the present results to contradict the proposed model. 

In fact, one can easily visualize the low-speed region, although it 
is cigar-shaped rather than hemispherical, plays a similar role to that 
of the hemispherical protuberance in the laminar flow. This is, however, 
in contrast to the description above in that the low-speed fluid is 
considered now to be the cause of the streamwise vortices rather than a 
consequence of them as viewed by the present author. Obviously, more 
work needs to be done before definite conclusions are drawn on this 
cause-and-ef feet relationship . 

IV. SUMMARY 

Detailed flow patterns associated with the bursting event in tur- 
bulent flows are constructed from the ensemble-averaged velocity field 
obtained by applying a conditional sampling technique to the data base 
generated by the large-eddy simulation of turbulent channel flow by 
Mom and Kim . 8 This is the first time such detailed flow patterns 
occurring during the bursting event have been obtained. In particular, 
during the ejection event, the streamlines and velocity vectors down- 
stream of the detection point exhibit clearly a pair of counterrotating 
streamwise vortices underneath the sweeping motion. Contours of 
streamwise vorticity in several streamwise locations suggest that the 
viscous-induced vorticity during the sweep event may be the source of 
the streamwise vortices. This is in contrast to the model suggested 
by Smith , 6 which conjectures that the streamwise vortices are formed 
by the low-speed region. Contours of the ensemble-averaged velocity 
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show that the low-speed fluid in the wall region is surrounded by the 
high-speed fluid, which is consistent with the description of Praturi 
and Brodkey . 1 3 
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FIGURE CAPTIONS 


FIG. 1 Streamlines in the (x-y) plane at z = 0. The upper figure 
is the conditionally averaged streamwise velocity at y + = 21. Note 
that each increment in x corresponds to Ax + - 62. 

FIG. 2 Streamlines in the (y-z) plane upstream of the detection point, 
x + ~ -124. Note that each increment in z corresponds to Az + - 15. 

FIG. 3 Streamlines in the (y-z) plane downstream of the detection 
point, x + ~ 124. 

FIG. 4 Velocity vectors in the (y-z) plane at x + - -124. 

FIG. 5 Velocity vectors in the (y-z) plane at x + z 124. 

FIG. 6 Contours of streamwise vorticity in the (y-z) plane at 

~ -124. 

FIG. 7 Contours of streamwise vorticity in the (y-z) plane at 

x + = 0. 

FIG. 8 Contours of streamwise vorticity in the (y-z) plane at 
x + = 124. 

FIG. 9 Contours of fluctuating streamwise velocity in the (x-y) plane 
at z = 0. 

FIG. 10 Contours of normal velocity in the (x-y) plane at z = 0. 
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FIG. 11 Contours of fluctuating streamwise velocity in the (x 
at y + ~ 12. 

FIG. 12 Contours of fluctuating streamwise velocity in the (y 
at x + ~ 124 


z) plane 


z) plane 
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Fig. 1 
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Fig. 2 
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Fig. 3 
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Fig. 4 
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Fig. 8 


22 





23 




/ S S' ~ Z 

////"<- N N 

/ / / / ' ^ " SZ -> N N \ 

X f ( / s -JZ ~N "N \ \ \ > 

J t / * / /✓ ^ — s -v 'N \ \ \ \ t 

^ / / s ' S s'' s 'v\NV\l N / / - 

/*— / ^////// / /«* *N \ NNU f / / //X 

' ^//// // ✓ /» ,- — ->\>||l#/////,/ 

« r >1 * i i / / // ////v^ggS 


</ ' /ZZZ^ / l N V. v ///<////'/// 

f « z//(iv v». "v-, *v •** / y S n **/ tut 



Fig. 10 
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Fig. 11 
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